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ABSTRACT 

Some of the processes Involved In the nucleatlon and growth of thin films 
were simulated by means of a digital computer. The simulation results were 
used to study the nucleatlon and growth kinetics resulting from the various 
processes. Kinetic results obtained for Impingement, surface migration. 
Impingement combined with surface migration, and Impingement combined with 
surface migration and with reevporatlon are presented. A substantial fraction 
of the clusters may form directly upon Impingement. Surface migration results 
In a decrease In cluster density, and reevaporation of atoms from the surface 
causes a further reduction In cluster density. 

1NTROOUC I ION 

The kinetics of nucleatlon and growth of thin films has been the subject 

of experimental and theoretical study for a long time. The standard 

theoretical treatment Is based on the rate equations derived by Zlnsmelster 1 
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and modified by others to take Into account the various processes 
Involved, e.g.. Impingement, surface migration, reevaporation, and cluster 
formation. These equations cannot, however, be considered a faithful 
description of the process, since, whereas they are deterministic In nature 
and yield a single solution, the various processes Involved In nucleatlon and 
growth are clearly random and cannot be expected to yield a repeatable result. 


1 


*NRC-NASA Research Associate. 


ORIGINAL. PAGE 19 
OF POOR QUALITY 

The study reported here Is an attempt of a different approach to the 

problem which takes Into account this random nature of the processes Involved, 

namely, computer simulation. Similar techniques have been recently applied to 
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the study of aggregation and diffusion-controlled cluster formation . In 
this paper only the general features of the results wir be presented. A more 
detailed account of the results will be given elsewhere. 

DESCRIPTION OF THE MODEL 

A square lattice was selected for the simulation, since by far most of 
the experimental studies of nucleatlon and growth used as a subtrate the (100) 
surface of NaCl, which Is of this type 10 ' 16 . The basic assumption made In 
the model are 

(1) Only single atoms Impinge on the substrate surface. 

(2) Only single atoms are mobile on the surface. 

(3) Mobile atoms migrate on the surface by Jumps to nearest-neighbor 
sites, and the jump frequency Is constant. 

(4) Only single atoms can reevaporate from the surface. 

(5) Only nearest- neighbor Interactions are taken Into account. 

(6) A cluster Is formed by two or more atoms occupying adjacent lattice 

sites. 

(7) No decomposition of clusters takes place. 

(8) The surface Is free from sites that cause preferential nucleatlon, 
e.g., defects or Impurities. 

(9) The simulation Is limited to a single monolayer, l.e., to the early 
stages of the Frank-van der Merwe and Stranskl-Krastanov mechanisms of film 
growth . 

(10) The boundary conditions taken are of an Infinite potential barrier, 
l.e., atoms cannot cross the boundary. 
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Most of these Assumptions were given and discussed by Logan 2 . No claim Is 
made as to the validity of any of them. (A study of the effect of at least 
some of them upon nucleatlon and growth behavior Is planned.) 

PROCEDURE 

The simulations were carried out on a minicomputer. A lattice of 400 
(20x20) points was the sample lattice. To check the effect of lattice size 
on the results, some test runs were also performed on a 1600-point (40x40) 
lattice. The results In these runs fell within the error range of the results 
obtained for the 400-point lattice. 

The x and y coordinates of the sites of atom Impingement were created 
by a random number generator with homogeneous distribution. A single atom on 
the surface can jump In one of the four directions shown In Fig. 1. The 
actual direction of each jump was also determined by means of a random number 
generator, which was programmed to yield as an output one of the numbers 1, 2, 
3, or 4. Once two atoms occupied nearest-neighbor positions on the lattice 
they were not allowed to change their positions. 

For each case studied eight simulation runs were performed, and the 
average as well as the standard deviation of the various parameters was 
evaluated. To check the sufficiency of the number of runs, 16 simulation runs 
were performed for some cases. The results were In the error range obtained 
with eight runs, and no reduction of the standard deviation resulted. 

RESULTS AND 0ISCUSS1ON 

To study the characteristics of each of the processes Involved In 
nucleatlon and growth, some of the simulations were first done separately, and 
then the combined process was studied. The simulations were used to estimate 
the number of atoms nucleated (l.e., atoms In clusters), the cluster density, 
and the size distributions of clusters. 


Impingement 

Some of the clusters are formed directly on Impingement without any 
surface migration. The physical process simulated In this case Is the one 
that takes place when a film Is deposited on a very cold substrate. The 
results for the number of atoms nucleated as a function of time are presented 
In Fig. 2. A substantial fraction of the Implnilng atoms are trapped In 
clusters, this fraction Increases as the process proceeds. This Is 
demonstrated better In Fig. 3 where the percentage of nucleated atoms out of 
the total number of Impinging atoms Is plotted versus time. Figure 4 shows 
the cluster density resulting from direct Impingement. The density first 
Increases, reaches a maximum at a coverage of *0.4, and then decreases due 
to the overlap of clusters. 

Surface Migration 

Next, the effect of surface migration following Impingement was studied. 
Physically, this would correspond to the case where a film Is deposited on a 
very cold substrate which Is later warmed up to allow migration of atoms on 
the surface. The results reported here (Figs ‘j and 6) are only for one value 
of coverage, e » 0.1. The process reaches saturation after 87124 jumps. 
Comparison of Fig. 6 with Fig. 4 allows one to estimate the fraction of 
clusters that are formed directly on Impingement relative to the total number 
of clusters obtained after migration Is completed. The number of clusters per 
site formed on Impingement at this coverage Is 0.013±0.002, the number of 
clusters per site following migration Is 0.03410.003. Thus, some 28±5 percent 
of the clusters are formed directly on Impingement, and the remaining 72±5 
percent are due to surface migration. 
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Combined Impingement and Migration 

This case of cabined Impingement and migration corresponds to the 

situation where a film Is deposited on a substrate that Is hot enough to allow 

surface migration yet Is too cold to allow reevaporation of adsorbed atoms. 

The simulation results, for the case where the Impingement rate Is 
-3 

2.5x10 atom/( site) (unit time) and the jump frequency Is 1 per unit time, 
are presented In Figs. 7 to 9. The effect of surface migration on growth 
kinetics Is best seen by comparing Fig. 3 with Fig. 8, where the percentage of 
nucleated atoms out of total number of Impinging atoms Is plotted versus time. 

The effect of surface migration on nucleatlon kinetics, l.e., on the 
cluster density. Is revealed by comparing Fig. 4 with Fig. 9: Surface 
migration causes a substantial decrease In cluster density. 

Reevaporation 

Finally, the effect of the reevaporation of adsorbed atoms from the 
surface on the nucleatlon and growth behavior was studied. Only one case Is 
presented here, namely, the one where the probability of each single atom 
reevapaoratlng Is 0.5. In other words, on each Jump every single atom has an 
equal chance of either landing on an adjacent lattice site or leaving the 
surface. 

The growth rate results for this case are given In Fig. 10. It Is worth 
noting that In this case almost all the atoms on the surface at any given time 
are bound In clusters, since the ones that are not have little chance of 
surviving there. This has the effect of about a fourfold slowing of the 
growth kinetics, as can be seen by comparing Fig. 10 with Fig. 7. The effect 
on nucleatlon kinetics Is shown In Fig. 11 where a threefold reduction In 
cluster density, compared with the case of no reevaporation (Fig. , can be 
seen. 
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CONCLUSIONS 


1 . Because of the random nature of the various processes Involved In 
nucleatlon and growth, a large scatter In results Is obtained. This Is 
reflected In the high values of the standard deviation. 

2. Depending on Impingement rate ard jump frequency, a substantial 
fraction of clusters may be formed directly upon Impingement. 

3. Surface migration results In a decrease In the cluster density. 

4. Reevaporation of atoms from the surface causes a further reduction In 
cluster density. 
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FIG. 3. Percent of atoms nucleated directly 
upon impingement Impingement rate, 

Z 5x10"* atom/(site x unit time). 



FIG. 4. Cluster density due to impingement Impingement 
rate, 2. 5x10”^ atom/(site x unit time). 
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FIG. 6. Cluster density due to surface miration. Jump frequency, 
1 per unit time* coverage, a 1. 
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FIG. 9. Cluster density due to combined impingement and 
surface migration. Impingement rate. Z 5x10“^ atom/ 
(site x unit time); jump frequency, 1 pr.r unit time. 
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FIG. 10 Nucleated atoms due to combined impingement, 
surface migration, and reevaporation. Impingement 
rate. 2. 5xl0” 3 atom /(site x unit time); jump fre- 
quency, 1 per unit time; coverage, a 5. 
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FIG. 11. Cluster density due to combined impingement, 
surface migration, and reevaporation. Impingement 
rate. 2. 5xl0 - ’ atom /(site x unit time); jump fre- 
quency, 1 per unit time; reevaporation probability, 


